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An XMM-Newton view of the 'bare' nucleus of Fairall 9* 
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ABSTRACT 

We present the spectral results from a 130 ks observation, obtained from the X- 
ray Multi-Mirror Mission-Newton {XMM-Newton) observatory, of the type I Seyfcrt 
galaxy Fairall 9. An X-ray hardness-ratio analysis of the light-curves, reveals a 'softer- 
whcn-brightcr' behaviour which is typical for radio-quiet type I Seyfert galaxies. More- 
over, we analyse the high spectral-resolution data of the reflection grating spectrometer 
and we did not find any significant evidence supporting the presence of warm-absorber 
in the low X-ray energy part of the source's spectrum. This means that the central 
nucleus of Fairall 9 is 'clean' and thus its X-ray spectral properties probe directly the 
physical conditions of the central engine. The overall X-ray spectrum in the 0.5-10 
keV energy-range, derived from the EPIC data, can be modelled by a relativistically 
blurred disc- reflection model. This spectral model yields for Fairall 9 an intermediate 
black-hole best-fit spin parameter of a — 0.39^ 
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1 INTRODUCTION 



Fairall 9 jFairallll 19771) is a type I Seyfer t galaxy at a redshift 
of z=0.047 (|Lauberts fc Valentiidll989h . It has a black- hole 
(BH) ma ss of (2.55±0.56) x 10 s M (based on reverberation 
mapping. IPeterson et alj2004h and an Eddington luminosity 
fraction Leoi/i/Edd of 0.05, assuming the bolometric lum i- 
nosity, LboI = 1.7 x 10 45 erg s"\ after I Woo fc Urrvl (|2002T ). 

F rom the early 70's Fa irall 9 (also known as ESO 113- 
IG 45. lHolmberg et allll978l ) has been observed severa l times 
in the X-ray band with Uhuru dForman et alJll978h, Ari el 
V dMcHardv et al.lll98ll). HEAO-1 A1 dWood et al.lll984) 
HEA O-1A2 jPiccinotti et all fl982l) Eins tein jPetre et all 
1 19841 ) and EXOSAT (|Morini et all Il986f ). The Advanced 
Satellite for Cosmology and Astrophysics (ASCA) revealed 
a broad iron line which was well mode lled by a relativis- 
tic line-emission model (|Revnolddfl997l ). Further observa- 
tions with the X-ray Multi-Mirror Mission-Newton (XMM- 
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Newton) observatory l|Gondoin et aljfeoOll ) unveiled the ex- 
istence of an absorption edge at 7.64 keV (in the source's 
rest-frame). In contrast to the ASCA observations, the iron 
line appeared to be narrow in the XMM-Newton data-set, 
suggesting that the reflection proc ess arises from material 
far fro m the putative central BH. iBrenneman fc Reynolds! 
(2009) reanalysed the same observation and found a third 
narrow emission-line at 6.78 keV (in the source's rest-frame) 
corresponding to ionized iron. 

Additionally, recent Suzaku obs ervations of Fairall 9 
showed a broad asymmetric iron line (|Schmoll et alJ feoOO). 
which is consistent with a disc-origin. A relativistically 
blurred disc reflection model, was fitted to the broad-band 
Suzaku spectrum (0.7-30 keV) of Fairall 9, yielding a con- 
straint for the BH spin parameter of a — 0.6 ± 0.07, where 
q = JM~ 2 (in the geometric unit system, iBardeen et al.l 
Il972h and —1 < a < 1, with J being th e angular momentum 
of the BH and M the BH mass. Finally. [Patrick et alj (|201ll ) 
reanalysed the same data using additionally data from the 
burst alert telescope (BAT) on board SWIFT satellite and 
derived a similar low value of a = 0.441q'j 4 . 

In the current paradigm, active galactic nuclei (AGN) 
are thought to consist of an accretion disc surrounding a 
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supermassive BH |Reeslll984l ). Consequently, most of the 
emitted energy that is released in the X-rays is the result of 
the conversion of gravitational potential energy from matter 
falling towards the BH. 

Thus, the X-ray spectroscopic properties of AGN can 
disclose valuable information about the physical properties 
of the processes operating in the immediate environment 
of the BH. In particular, broadening of the iron emission- 
lines, formed in the inner part of the accretion disc, can pro - 
vide constraints about the BH spin parameter (|Millerll2007h . 
Since the X-ray nucleus in Fairall9 is 'clean', la cking X-ray 
warm-absorption features (|Gondoin et al-lfeOQll ). it offers an 
ideal environment to investigate the physical properties of 
the central region. 

However, there are alternative scenarios that can ex- 
plain the broadening of the iron lines. Complex absorption 
processes from warm-absorbers can cause an apparently sim- 
ilar iron line p rofile to the profiles caused by strong gravity 
environments l|Reeves et al.ll2004l ; lTurner et al.ll2005l ). More- 
over, it is reasonable to expect several weak helium- and 
hydrogen-like emission-lines from distant material around 
the same energy as the iron line (i.e. 6.4-6.7 keV energy- 
range in the sourc e's rest-frame) , caus ing si milar apparent 
broad ening effects (|Nandra et al.| [2007). Also lNandra et all 
l|2007l ) showed that intermediate ionization species of the 
iron itself, can produce asymmetric broad line-profiles ex- 
tending up to 6.97 keV (in the source's rest-frame). Fi- 
nally, blueshifted or redshifted spectral-features from ejected 
or infalling materiel, close to the nucleus, can appear in 
the energies around the iron line, br oadening its profile 
i|Yaqoob et al.lll999l ; [Turner et al.ll2002l ). 

In this paper we report the analysis results of an 
XMM-Newton observation of Fairall 9 performed in Decem- 
ber 2009. Initially, in Sect. [2] we present the data-reduction 
procedures for the X-ray data from the European photon 
imaging camera (EPIC), consisting of the pn-charge coupled 
device (pn-CCD) and the two metal oxide semi-conductor 
(MOS)-CCDs. In the same section, we present the data- 
reduction details for the soft X-ray data from the reflecting 
grating spectrometer (RGS). Then, in Sect. [3] we show the 
EPIC light-curve products together with a hardness-ratio 
analysis. In Sect. [4] we perform a spectral analysis of the 
RGS data and in Sect. [5] we investigate the spectral proper- 
ties of the EPIC data. A discussion together with a summary 
of our results can be found in Sect. [6] The cosmological pa- 
rameters used throughout this paper are: Ho = 70 km s _1 
Mpc -1 , Oa=0.73 and f2 m =0.27, yielding a luminosity dis- 
tance to Fairall 9 of 208.5 Mpc. 



2 OBSERVATIONS AND DATA-REDUCTION 

2.1 EPIC-MOS and EPIC-pn data-reduction 

Fairall 9 has been observed by XMM-Newton (Obs-ID: 
0605800401) from 2009 December 9, 19:56:35 (UTC), to 
2009 December 11, 08:04:07 (UTC) (on-time: 130052 s). The 
pn camera was operated in Prime Small Window mode and 
the two MOS cameras in Prime Partial W2 mode. Medium- 
thickness aluminised optical blocking filters were used for 
all EPIC cameras to reduce the contamination of the X-rays 
from infrared, visible, and ultra-violet light. 



The EPIC raw-data are reduced wi th the XMM-Newton 
Scientific Analysis System (SAS) (|Gabriel et alj|2004l ) 
version 10.0.2. After reprocessing the pn and the two MOS 
data-sets with the epchain and the emchain SAS-tools re- 
spectively, we perform a thorough check for pile-up using the 
task epatplot. The pn data do not suffer from pile-up, thus 
the source and the background count-rates are extracted 
from a circular region of radius 31.5 arcsec and 35 arcsec, 
respectively. Conversely, both MOS cameras appear to be 
slightly piled-up. Therefore, to be conservative, the final 
source count-rates are extracted from annuli of the same 
outer radius 35 arcsec and inner radius 2.75 arcsec and 2.5 
arcsec for the MOS 1 and the MOS 2, respectively. The back- 
ground count-rates are extracted from two circular regions of 
3.33 arcmin and 3.42 arcmin for the MOS 1 and the MOS 2, 
respectively. Note that these apertures yield the optimum 
signal-to-noise ratio above 7 keV. We verify that the result- 
ing light-curves are not affected by pile-up problems. 

For the production of the light-curves at a given energy- 
range, we select events that are detected up to quadru- 
ple pixel-pattern on the CCDs i.e. PATTERN<12. The 
corrected background-subtracted light-curves of the source 
are produced using the SAS-tool epiclccorr. Based on the 
pn background light-curve, in the 0.5-10 keV energy-range 
(Fig. [U grey-points corresponding to the right-axis), an in- 
creased background activity was registered during the first 

4 ks and the last 7.9 ks of the observations and thus, data 
from those periods are disregarded from our analysi^J. This 
leaves us with 118152 s clean exposure-time. 

For the production of the spectra we use the SAS-task 
evselect and we select for the pn camera the pulse-invariant 
(PI) channels from to 20479, having a spectral-bin-size of 

5 eV. For the two MOS cameras we select the PI channels 
from to 11999 with a spectral-bin-size of 15 eV. For the 
spectral analysis we have set more stringent filtering criteria 
from those in the case of the light-curves. For the pn data 
we allow events that are detected up to double pixel-pattern 
on the CCD i.e. PATTERN <4 and we exclude all the events 
that are at the edge of a CCD and at the edge to a bad pixel 
i.e. FLAG=0. For the MOS data we allow, as in the case 
of the light-curves, PATTERN<12 but we restrict to those 
events having also FLAG=0. After filtering the events, we 
calculate the area of the source and the background annular 
regions using the task backscale and finally we compute the 
detectors' response matrices and the effective areas using 
rmfgen and arfgen, respectively. 

The X-ray spectral fitting analysis (Sect. [5t has been 
performed by the high-level task XSPEC (|Arnaudl 19961 ) ver- 
sion 12.6.0, of the HEAsoft version 6.10. We group the PI 
channels of the so urce-spectra, mak ing use of the task grppha 
of the FTOOLS l|Blackburnlll995l ) version 6.10. The MOS 
source-spectra are grouped in bins of four (A_E=60 eV) and 
eight (AiJ=120 eV) spectral-channels for the energy-ranges 
of 0.5-8 keV and 8-10 keV, respectively. The pn source- 
spectra are grouped in bins of four (AE =20 eV) spectral- 



1 Despite the fact that the 'high' background-states are still min- 
imal, in terms of count-rate, with respect to the source count-rate 
in 0.5—10 keV (Fig. [pj black-points corresponding to the left-axis), 
they consist of events having energies more than 2 keV. At these 
energies the source count-rate is of the order of 1.5 counts-s -1 . 
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channels for the energy-range of 0.5-10 keV. All the grouped 
spectral-bins contain at least 20 photons. 



2.2 RGS data-reduction 

The RGS raw-data are processed following the standard 
data-reduction threads of SAS. We use the tool rgsproc to 
extract calibrated first-order spectra and responses for the 
RGS 1 and RGS 2 cameras. In general, XMM-Newton RGS 
observations can be affected by high particle background 
periods during parts of the XMM-Newton orbits, mostly 
caused by Solar activity. The high energy-band is the one 
being mostly affected by background flares, and high energy 
RGS photons are dispersed over the CCD-9 chip. We there- 
fore extract the background light-curve of the CCD-9 chip 
and select, as good-time- intervals of the processed obser- 
vation, only those during which the background count-rate 
deviates by less than two standard deviations from the av- 
erage background count-rate of each observation. The final 
(i.e. after cleaning for high-background time-intervals) ex- 
posures of the RGS 1 and RGS 2 spectra are 129.6 ks and 
129.8 ks, respectively. 

Both RGS 1 and RGS 2 spectra are grouped at a res- 
olution of 8 mA (~0.16 eV at 0.5 keV), so allowing 8 
spectral- bins for each RGS resolution element. For both the 
RGS 1 and the RGS 2, the photon count-rate of the back- 
ground spectrum becomes comparable to the background- 
subtracted source count-rate, at energies less than 0.35 keV 
(A > 35 A), and greater than 1.8 keV (A <7 A). Moreover, 
due to failures of two different read-out detector chips, early 
in the life of the mission, both the RGS 1 and the RGS 2 
lack response in two different spectral intervals of 0.9-1.2 
keV and 0.5-0.7 keV, respectively. We therefore consider the 
following two spectral intervals for spectral fitting purposes: 
[(0.35 - 0.9) U(1.2 - 1.8)] keV and [(0.35 - 0.5) (J(0.7 ~ 1-8)] 
keV, for RGS 1 and RGS 2, respectively. 

For the spectral fits the Sherpa modelling and fitting ap- 
plication of the Chandra Interactive Analysis of Ob- 
servation (CIAO) software version 4.2, has been used to 
fit simultaneously the RGS 1 and RGS 2 spectra of Fairall 9. 
In addition, the spectral fitting analysis results have been 
cross-checked and verified with XSPEC, but the latter val- 
ues are less conservative and thus we report only the results 
coming from Sherpa. 



3 X-RAY LIGHT-CURVES 

The combined X-ray light-curves from the EPIC-pn and the 
two EPIC-MOS cameras in 0.5-1.5 keV, 2.5-4 keV, and 5-10 
keV energy-ranges, are shown in Fig. [21 A steady count-rate 
increase can be noticed in all three energy-bands, of the or- 
der of 30 per cent, 20 per cent, and 10 per cent, respectively. 
Additionally, the fractio nal variability amplitu des (corrected 
for the photon noise e.g. IVaughan et al. 2003), for the three 
bands are 7.3±0.1 per cent, 5.7±0.6 per cent, and 2.2±2.2 
per cent, respectively. Both the overall count-rate and the 
fractional variability amplitude increase with decreasing en- 
ergies. 
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Figure 1. The Fairall 9 combined light-curve (black-points, cor- 
responding to the left-axis) and the background pn light-curve 
(grey-points, corresponding to the right-axis) in the 0.5-10 keV 
energy-band. The background shows increased activity during the 
first 4 ks and the last 7.9 ks. 
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Figure 2. The combined EPIC-MOS 1,-MOS 2,-pn X-ray light- 
curves of Fairall9 in the 0.5-1.5 keV, 2.5-4 keV, and 5-10 keV 
energy-bands, in bins of 100 s. Note the difference in the count- 
rate scale on the vertical axis. 
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that we get from the data-set is 18.71, the probability of 
getting such a value by chance alone is only 1.42 x 10 -14 . 

For HR2 the best-fit model gives a slope of ( — 1.5 ± 
0.5) x 1CT 2 having a 95 per cent c.i. of (-2.6 x 1(T 2 , -0.5 x 
10 -2 ). The value of f-statistic for this data-set is 11.95 yield- 
ing a chance coincidence probability of only 7.86 x 10~ . 
Since in both cases the zero value lies outside the estimated 
c.i., the null-hypothesis (i.e. the slope is equal to zero) is 
rejected at 5 per cent significance level. 

This result implies that Fairall 9 shows significant X-ray 
spectral variations, becoming 'softer-when-brighter'. This 
beh avior is considered typical for ra dio-quiet type I Seyfert 
(e.g. ISobolewska fc Papadakis!l2009l ). 



4 THE RGS SPECTRUM OF FAIRALL 9 

We search the Fairall 9 RGS 1 and RGS 2 spectra (in their 
available wavelength intervals) for narrow absorption and 
emission features. We inspect locally 2 A wide spectral re- 
gions, and we model the local continua with simple power- 
laws attenuated, in this line of sight, by the Galactic neu- 
tral absorption Nn = 3.2 x 10 20 cm~ 2 (estima ted using the 
FTOOL nH, after iDickev fc Lockmanl Il990l 'l. When both 
the RGS 1 and the RGS 2 are available we fit simultane- 
ously the same 2 A regions from the two spectrometers. In 
Fig. [5] we show the 18-24 A portion of the RGS 1 spectrum 
of Fairall 9, together with its best-fit continuum model (red- 
line) and the corresponding ratio plot (data/model). Note 
that in this section all the centroid-energies of the emission- 
absorption-lines refer always to the observer's frame. 

From the ratio plot of Fig. [S] we find no significant 
evidence for absorption-lines imprinted either by intrinsic 
(i.e. nuclear) or by intervening (i.e. interstellar medium 
(ISM) or intergalactic medium) material down to 4 stan- 
dard deviations (below the fluctuations of the continuum) 
line equivalent- width (eqw) sensitivity of 45 mA. The only 



exception is a relatively strong local O I Ka line (eqw= 36_jj 
mA), imprinted by our Galaxy ISM at E=0.5273±£ keV 
(23.51^0 06 A), and detected at a confidence level of 2.9 stan- 
dard deviations (below the fluctuations of the continuum). 
Assuming a relative fraction (with respect to the ensemble 
of the oxygen species) of O I of the order of unity, i.e. neu- 
tral material, this gives an average ISM metallicity along 



The value of t-statistic is t2i, 0.025 = 2.08 
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3.1 Hardness-ratio analysis 

Observing different variability amplitudes in different 

energ y-bands is usually in dicative of spectral variations 

(e.g. iBrinkmann et al.ll2007h . After binning the light-curves 

in bins of 5 ks, we estimate the hardness-ratios (5-10 

keV)/(0.5-1.5 keV) (HR1) and (5-10 keV)/ (2.5-4 keV) 

(HR2) versus the overall count-rate in 0.5-10 keV (Fig. [3]). 

We see a decreasing trend implying that the source X-ray 

spectrum becomes softer when the source gets brighter. 

In order to check the significance of such a trend, we fit 

to the 23 hardness-ratio points a linear model and we check 

how significantly the value of the slope differs from that of 

zero. For HR1 the best-fit-model gives a slope of (—3.3 ± 

0.5) x 10~ 3 with a 95 per cent confidence interval (c.i.) of 

(-4.3 x 10~ 3 ,-2.3 x 10~ 3 j!|. Since the value of i-statistic > °- 66 
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Figure 3. Hardness-ratios in (5-10 keV)/(0.5-1.5 keV) and (5- 
10 keV)/(2.5-4 keV) versus the total count-rate in 0.5-10 keV, 
binned in 5 ks time-bins. The solid-line represents the best-fit 
linear model. The negative slope implies a 'softer-when-brighter' 
X-ray spectral behavior. 



the line of sight of the order of 20 per cent that of the So- 
lar value. Moreover, we find evidence of two additional local 
absorption-lines OvnK/3 (eqw= 39t*o m ^) an d OvuKa 
(eqw= 38±^ mA) at E=0.6640t°;°°'j' 3 keV (18.67l°;°4 A ) 
and E=0.5746± o ;^^ keV (21.58l^ 2 A), respectively. These 
lines are detected at a confidence level of 3.8 and 3.2 stan- 
dard deviations (below the fluctuations of the continuum) 
and they are thought to be created by an ubiquitous lo- 
cal absorber (seen along virtually all lines of sight), which 
could be hot gas in an extended corona of our galaxy, or a lo- 
cal warm-hot intergalactic medium filament permeating our 
own L ocal Group and/or Milky Way, or both (|Nicastro et al.l 
2002). However, from the ratio plot of Fig. [5] (dashed-lines 
at 18.67 A and 21.58 A) we notice that each of these two 
absorption-lines coincides with effective area features, i.e. 
instrumental features. 

Additionally, we find evidence for the presence of few 
narrow (i.e. unresolved) emission-lines close to the systemic 
redshift of Fairall 9. In particular, we identify two of the 
three lines of the OvuKa emission triplet, the resonant 
(r) and the forbidden (f) lines (Fig. [4] solid-lines). Dur- 
ing the fit, the three lines (r,f and intercombination (i)) 
have positions linked to their rest-frame line energy-ratio 
and only the energy of O VII Ka-r is left to vary freely. The 
rest-frame wavelengths of these lines are 21.602 A (r), 21.8 
A (f), and 22.1 A (i). Unfortunately, the OvuKa-i falls 
into an instrumental feature of the RGS 1 (the RGS 2 is 
blind at these wavelengths) and thus no density diagnos- 
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tics are possible for the X-ray emitting gas. The VII Ka-f 
line (eqw= 180±|| mA) at E=0.5379 keV (23.05 A) and 
O vnKa-r line (eqw=42l^ m A) at E=0.5502± 0.0004 keV 
(22.53±0.02 A) are detected at a confidence level of 4.5 and 
1.9 standard deviations (above the fluctuations of the con- 
tinuum) at a redshift of z=0.043±0.001, implying an outflow 
velocity relative to the optical systemic redshift of Fairall 9 
of -1200±300 km s' 1 . The OviiKq (f/r) 4.3±|;| intensity 
ratio strongl y suggests the photo-ion ization nature of the 
emitting gas l|Porquet fc Dubaull2000h . 

Despite the absence of individually significant 
absorption-lines, we check further for the presence of 
any ionized gas in the source. To this end, we jointly fit 
the RGS 1 and RGS 2 spectra with a model consisting of a 
power-law plus a warm-absorber attenuated, in this line of 
sight, by the Galactic neutral absorption Nh- We use the 
self-consistent c ode Photoionized A bsorption Spectral 
Engine (PHASE) (|Krongold et alj|2003l ), that calculates the 
absorption caused by ionized plasma. The code has four 
parameters: the equivalent hydrogen column density of the 
absorber, the ionization parameter U of the gas (defined as 
the ratio between the density of ionizing photons and the 
electron density of the gas), the turbulent velocity, and the 
redshift of the absorber. 

We freeze the absorber's redshift to the systemic red- 
shift of Fairall 9, and its turbulent velocity to v tur b = 100 km 
s _1 , while we leave the equivalent hydrogen column density 
and the ionization parameter U of the absorber free to vary 
during the fit. The best-fit value that we obtain for log 10 (U) 
is — 1-llto'g is rather low for typical type I Seyfert warm- 
absorbers. Furthermore, although we can obtain just an up- 
per limit for the equivalent hydrogen column-density of the 
absorber, this limit is rather tight: log 10 (JVh cm -2 ) ^ 20.36 
(the errors correspond to one standard deviation for the two 
interesting parameters). In Fig. [S]we show the 68.27, 95.45 
and 99.73 per cent confidence contours for these two param- 
eters. Clearly, a significant amount of intrinsic absorbing 
material along our line of sight to Fairall 9 is not required 
by the RGS data. 



5 EPIC X-RAY SPECTRAL ANALYSIS 

For the X-ray spectral results, all the errors on the best-fit 
parameters indicate their 90 per cent confidence range, cor- 
responding to a A\ 2 of 2.706, unless otherwise stated. Dur- 
ing the fitting procedures, we keep all the model-component 
parameters between the two MOS and the PN spectra tied 
together except from the normalizations (for the two MOS 
spectra they are still tied together). The model-flux is modi- 
fied using the photoelectric absorption model wabs having a 
fixed interstellar column-density of Nh ■ Note, that the best- 
fit values for the centroid energy of the emission-lines are 
given in the source's rest-frame and the horizontal axis in 
the spectral-plots always refers to the energies in the ob- 
server's frame. 



5.1 A first look at the broad-band spectrum 
(Model 1) 

In order to get a phenomenological view of the overall spec- 
tra we initially fit a redshifted power-law (Model 1) to the 



0.015 



a 

3 

o 
U 




0.005 



19 



20 21 22 
Wavelength (A) 



23 



24 



Figure 4. Spectral fit in the 18-24 A wavelength-band of the un- 
binned (i.e. AA = 8 mA) RGS 1 spectrum. The spectral model is a 
power-law attenuated Galactic neutral absorption (red-line). No- 
tice that the deep absorption-like features, appearing in this spec- 
tral interval, are instrumental features of the RGS 1 dispersion 
CCD detectors (i.e. dead pixel/columns), and are modelled by the 
RGS 1 effective area. The only real absorptions features are the 
weak absorptions-lines (dashed-lines) created locally (i.e. Milky 
Way and/or in the Local Group) OlKa, OvnK/3, and OvoKa. 
The only real emission features present are the O VII Kct-r and 
O VII Koj-/ emission-lines (solid-lines) detected at the systematic 
redshift of Fairall 9. 



data in the 3-5 keV and 7-10 keV energy-bands. The \ 2 
value of the best-fit model is equal to 419.41 for 388 d.o.f. 
yielding a best-fit photon spectral index of F = 1.66 ± 0.02. 
Figure [6] shows this best-fit model extrapolated over the 
energy-range of 0.5-10 keV together with the respective ra- 
tio plot. 

From the ratio-plot of Fig. [5] we can readily distinguish 
two very distinct spectral features. Below 2 keV, we notice 
a strong flux-excess on the top of the best-fit power-law ex- 
trapolation. Then, in the 5.7-8 keV energy-range a strong 
emission-line appears around 6.1 keV followed by an excess 
around 6.7 keV. In addition to those features, significant 
broad residuals appear above 7 keV as well as a feature at 
6.8 keV suggestive of an absorption edge. 



5.2 Spectral fitting in the 3—10 keV energy-range 

Since Model 1 (Fig. [6} gives an adequate description of the 
data in the 3-5 keV and 7-10 keV energy-ranges, we keep it 
as a baseline model, for the following spectral-analysis over 
the 3-10 keV energy-range. The best-fit parameters for the 
following X-ray spectral models are given in Table [1] 
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logio(U) 

Figure 5. The 68.27 per cent, 95.45 per cent, and 99.73 per 
cent confidence contours for the best-fit parameters, log 10 (U) and 
logig (^H cm_ 2 ) indicated by the cross, for the best-fit power- 
law plus warm-absorber model to the RGS 1 and RGS 2 spectra. 
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Figure 6. Model 1: Redshifted power-law fit in the 3—5 keV and 
7-10 keV energy-ranges, indicated by the grey-areas. The best-fit 
model is extrapolated over the energy-range of 0.5-10 keV 



model to account for the associated Fe K/3 lin^E For the lat- 
ter, we fixed the centroid energy at 7.06 keV (in the source's 
rest-frame) and its norma lization at a value eq ual to 0.16 of 
the Fe Ka normalization (|Molendi et alj|2003r i. 

We model the absorption-edge and the associated 
reflection-continuum as the reflection of an exponentially 
cut-off power-law spectr um from a neutral disc using the 
XSPEC model pexrav (|Magdziarz fe Zdziarskil 1 19951 ). In 
order to derive the reflected continuum only, we restrict 
the range of the relative-reflection scaling-factor to nega- 
tive values. The energy cut-off is fixed at £^=200 keV (in 
the source's rest-frame), the abundance of elements heavier 
than helium is set equal to that of the Solar value, the iron 
abun dance is set equal 0.8 of the Solar value |Schmoll et al.l 
2009), and the inclination angle is fixed at <f> = 40 degrees. 
Note that we do not tie together the photon index of the 
pexrav model with that of the redshifted power-law. If the 
reflecting material is associated with the distant molecular 
torus, which is expected in Seyfert galaxies, it will proba- 
bly respond to the average X-ray source-spectrum, whose 
slope may be different than the spectral slope during this 
XMM-Newton observation. 

The best-fit model (Model 2) is plotted in the top panel 
of Fig. [71 having a x 2 value of 600.29 for 538 d.o.f. The best- 
fit centroid energy of the Fe Ka line is 6.42T0.01 keV. The 
best-fit values for the relative reflection and the power-law 
photon index, from the pexrav model, are 7?pex = 0.87^ 'o6 
and Fpex = 1.87 ± 0.03, respectively. Finally, the best-fit 
power-law photon index, coming from the redshifted power- 
law model-component is T = 1.85* 'o2- From the ratio plot 
of the top panel of Fig. we see that even after including 
the two narrow Fe lines and the reflection component, three 
broad emission components still remain around 5.9 keV, 6.1- 
6.7 keV, and 6.9-7.7 keV (in the observer's frame). 

Note that the photon indices coming from the pexrav 
and the power-law model-components are consistent within 
their errors. Tying the two photon-indices together and re- 
peating the fit yields a value of V 
tent with the previously derived values. 



1.871q'q2> being consis- 



5.2.2 Adding a relativistic disc-line (Model 3) 

In order to quantify further the aforementioned spectral- 
excesses, we fix the parameters of Model 2 to their best- fitted 
values and we add the relativistic disc-line model of iLaorl 
|l99ll ). denoted as laor in XSPEC, (Model 3). The latter 
model-component describes the line-profile emitted from an 
accretion disc around a rotating black hole including the 
strong relativistic effects that are produced in the innermost 
regions of the accretion disc. 

During the fit, we fix the parameters of the inner and the 
outer radius of the accretion disc to 6 and 300 gravitational 



5.2.1 Iron line and reflection (Model 2) 

To account for the emission-line at 6.1 keV (in the observer's 
frame), we consider a redshif ted zero-width Gaus sian- 
function model. Foll owing both iGondoin et all (|200lh and 
ISchmoll et all (|2009l ) and due to the fact that this line ap- 
pears to correspond to the Fe Ka line from neutral iron, hav- 
ing a centroid energy of 6.4 keV in the source's rest-frame, 
we add further a redshifted zero-width Gaussian-function 



3 The exact energy position of these two iron lines depends on the 
ionization state of the Fe element in th e source. Assuming neu- 
tral atoms (after Kaastra & Mcwc 1993) the Fe Ka is a doublet 
composed of a 6.391 keV and a 6.405 keV line, as well as the Fe 
K/3 line consisting of a 7.057 keV and a 7.058 keV line. Since the 
energy resolution of the EPIC cameras at 6-7 keV energy-range 
is around 150 eV (full width half maximum), it is impossible to 
resolve them and thus we fit them with a zero-width redshifted 
Gaussian- function. 
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radii, r g , respectively. From the bottom panel of Fig. [7] we 
can see that Model 3 gives a very good description of the 
excesses in the 6.1-7.7 keV energy-range (in the observer's 
frame). The \ 2 value of this fit is 559.18 for 536 d.o.f. The 
best-fit value for the centroid energy of the relativistic line 
is 6.75±0.18 keV. Also, the best-fit values for the emissivity 
index and the inclination angle are q = 2.46lg'37 and 6 = 
84±i 7 de grees, respectively. Finally, the power-law photon 
index is steeper than that of Model 2 and equal to T = 
1.97^Q Qg. In Fig.[8]we show the spectral components of both 
the best-fit Model 2 (dashed-line) (Sect. l5.2.Tj) and the best- 
fit Model 3 (solid-line) to the pn data, depicting the relative 
flux contribution of each component to the X-ray spectrum. 

For this model, we can now estimate the eqw of the 
narrow and broad emission-lines, having as an underlying- 
continuum the contribution from all other components of 
Model 3. The eqws for the Fe Kq and Fe K/3 lines that we 
get are 103±19 eV and 19±4 eV, respectively (mean values 
from the pn and the MOS data). The eqw of the relativistic 
disc-line is 362lgjg eV and 1621^2 e V, for the pn and the 
MOS data, respectively. 

Theoretically, it is expected that the eqw of the narrow 
emission-lines should be tightly connecte d with the amoun t 
of the continuum reflection component l|Matt et al.lll99ll ). 
In fact, one can consider a linear correlation between the 
two parameters with the reflection strength, R = (Q. 
is solid angle that the reflector subtends) bei ng equal to 
unity when the Fe Kq lin e's eqw is 120 eV (e.g. iMatt et al.l 
Il99ll ; iBianchi et ai] 120091 ). For the best-fit eqw- value of Fe 
Ka line, that we have found previously, the expected re- 
flection should be R = 0.86 ± 0.17. This value is entirely 
consistent with the reflection fraction estimates mentioned 
above (Model 2). Based on the estimated reflection value, 
the average solid angle that the distant reflector should sub- 
tend is (1.72±0.34)tt sr. 

Note that if we replace the relativistic disc-line model 
with a narrow emission-li ne at 6.78 keV (in the source 's 
rest-frame), as observed bv lBrenneman fe Revnoidi (2009), 
the broad excess between 6.1-6.7 keV (in the observer's 
frame) can not be modelled. Instead, we need at least three 
emission-lines to account for this excess resulting to an over- 
parametrised physical model. 



5.2.3 Ionized disc reflection with relativistic blurring 
(Model 4 ) 

The simple relativistic disc-line component of Model 3 seems 
to give a good morphological description about the nature of 
the spectral excesses of Model 2. Nevertheless, the inclina- 
tion angle, derived by the laor model, seems to be unphys- 
ically high for a type I Seyfert galaxy. Also, the centroid en- 
ergy of the relativistic line strongly indicates that the inner 
disk is moderately ionized, yet no ionized Compton reflec- 
tion component is included in Model 3. In this section we 
consider a relativistic self-consistent spectral-model, using 
as an input the reflection coming from a constant density il- 
luminated disc. This model (Model 4) cons i sts of the X-ray 
ionized-reflection model of iRoss fe Fabian] l|2005h . denoted 
as ref lionx in XSPEC, con volved with the relativist i c blur - 
ring multiplicative model of lBrenneman fe Reynolds! l|2006h . 
denoted as kerrconv in XSPEC. 

The ionized-reflection model takes into account several 
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Figure 7. Spectral fits in the 3—10 keV energy-range. The vertical 
labels indicate the energy in the source's rest-frame. (Top panel) 
Model 2: Two narrow Fe lines -Fe Ka (free), Fe K/3 (fixed)- 
, a redshifted power-law, and a reflection component. The in- 
set shows the model together with the actual spectral channels 
of the pn and MOS detectors, separated by 10 eV and 15 eV, 
respectively. (Lower panel) Model 3: Two narrow Fe lines — Fe 
Ka (fixed), Fe K/3 (fixed) — , a redshifted power-law, a reflection 
component (fixed), and a relativistic disc-line. 



spectral signatures of ionized reflection in the source's rest- 
frame including: the Compton-hump at 30 keV, a strong 
ionized iron line at 6.7 keV, multiple emission-lines on top of 
a soft continuum between 0.3-3 keV, and a further extreme- 
UV emission hump peaking at about 60 eV. The blurring 
model smears out the spectrum due to the strong gravity 
effects in a fully relativistic way and it allows the BH spin 
parameter, a, to be a free fitting parameter. 

During the fitting procedure, we keep all the param- 
eters of Model 3 (i.e. the Fe Ka and the Fe K/3 lines, as 
well as the pexrav model-component) fixed at their best- 
fit values and we replace the disc-line component with the 
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Figure 8. Spectral components for the best-fit Model 2 (black 
dashed-lines) and best-fit Model 3 (black solid-lines) to the pn 
data. The red lines correspond to the overall best-fit theoretical 
models respectively. 



kerrconvxref lionx. We assume a disc of uniform emissiv- 
ity extending from the radius of the innermost stable cir- 
cular orbit, -Risco, out to 400 -Risco- The photon index 
of the ionized-reflection model-component, F, is tied to the 
photon index of the underlying redshifted-power-law model- 
component, representing the continuum. 

The best-fit Model 4 gives a statistically acceptable fit 
for the X-ray EPIC-data of Fairall 9 (Fig. EJ| yielding a \ 2 
value of 578.88 for 535 d.o.f. The best-fit model parameters 



0.20 



are: power-law photon index of Y = 1.84 



+0.02 



emissivity in- 



dex of q = 2.07^q ionization parameter of £ = 18. SJ^ 1 ^ 1 
erg cm s _1 , inclination of 6 — 50 jl^f degrees, and BH spin- 
parameter of a = 0.43lo Qg. 

Note that our data do not allow us to discern between 
uniform and non-uniform emissivity profiles. If we consider 
a kerrconv model-component with a non-uniform emissivity 
profile, we get two emissivity indices of qi = 2.65l 2< ?9 and 



q 2 = 1.58±g;g|, and a break-radius of r br = 20+^' r g . The 
two indices are consistent with the value derived previously 
for the case of uniform emissivity. 



5.3 The broad-band spectrum, 0.5—10 keV: 
Ionized disc-reflection with relativistic 
blurring (Model 4 (broad-band)) 

As we have mentioned above, at soft energies (i.e. below 2 
keV) in addition to the power-law continuum, the ionized- 
reflection model predicts the presence of multiple emission- 
lines that are expected to be almost entirely blurred, due to 
the strong relativistic effects that operate in the vicinity of 
the central BH. In fact, it has been suggested that the so- 
called 'soft-excess' emission in AGN could be the result of X- 
ray reflection from th e inner parts of the accretion disc (e.g. 
ICrummv e t al. 2006). For this reason, we examine whether 
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Figure 9. Spectral fits in the 3-10 keV energy-range of Model 
4. The vertical labels on the vertical lines indicate the energy in 
the source's rest-frame. Two narrow Fe lines -Fe Ka (fixed), Fe 
K/9 (fixed)-, a reflection component (fixed), a redshifted power- 
law, and an ionized disc-reflection convolved with a relativistic 
blurring model. 



Model 4 is able to reproduce the broad-band X-ray spectrum 
of the source in the 0.5-10 keV energy-band. 

To this end, we extend the energy-range down to 0.5 
keV and we fit again Model 4 to the pn and MOS spectra 
simultaneously, having again the Fe Ka and the Fe K/3 lines, 
as well as the pexrav model-component fixed at their best- 
fit Model 3 values. The \ 2 value of the fit is 1052.63 for 749 
d.o.f. and the best-fit-model is shown in Fig. 1101 The best- 
fit model parameter values are: a power-law photon index 
of T = 2.01±n an emissivity index of q = l.98±° ° 9 5 , an 

-tie er S cm s_1 



ionization parameter of £ = 15. 82^3'^ erg cm s _1 , and an 
inclination angle of 6 = 64^g degrees. Finally, the best-fit 
BH spin parameter is 01 — 0.391q'3q. We outline the fact 
that the value of a is very sensitive to the rest of the fit- 
ting parameters. Even small deviations of the rest of the 
fitting parameters, around their best-fit values and within 
their quoted errors, can change significantly the result of a. 

The value of the x 2 fit-statistic for this best-fit model is 
rather high. Strictly speaking, for the given d.o.f., the null 
hypothesis probability is extremely small, 1.3xl0 -12 , im- 
plying that the particular best-fit model is highly unlikely 
to be a good representation of the underlying spectrum in 
the 0.5-10 keV energy-range. However, we believe that there 
are good reasons that this may not be the case due to the 
larger collecting area at lower energies. The quality of the fit 
is mainly determined by the data below 2 keV which have 
higher signal-to-noise ratio, than those above 2 keV, and 
thus bigger statistical- weight, during the estimation of \ 2 ■ 
At energies below 2 keV the cross-calibration uncertainties 
between the pn and the MOS detectors, as well as the indi- 
vidual calibration uncertainties of eac h detector separatel y 
are of the order of 8-10 per cent (e.g. Ma teos et aL| [2009). 
The dashed- and the dotted-lines in the ratio plot of Fig. [10] 
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indicate the 8 and 10 per cent residual-amplitude areas. We 
see that all of the ratio estimates, below 2 keV, are within 
these areas. What is even more interesting is that the pn and 
the MOS residuals do not typically agree with each other, as 
we would expect if indeed the model did not provide an ac- 
ceptable representation of the observed spectra, suggesting 
that calibration inaccuracies dominate the x 2 fit-statistic. 
Finally, in order to limit the systematic calibration uncer- 
tainties to those originating only from the pn instrument, 
we ignore the MOS data and we fit the pn spectra alone. 
In this case, the \ 2 °f the overall fit is better (572.91 for 
469 d.o.f.) but now due to the lower signal-to-noise ratio, 
since the MOS data have been excluded, we can not con- 
strain adequately the fitting parameters e.g. the black-hole 
spin parameter is a = 0.99^ 98 and the inclination angle 
9 = 90Z 81 degrees. 

Furthermore, as we have shown in Sect. 13.11 Fairall9 
does show flux variations during this XMM-Newton obser- 
vation which are associated with spectral variations. In fact, 
the largest amplitude variations appear in the soft-band i.e. 
0.5-1.5 keV. Therefore, it is not surprising that the fit- 
quality of (any) physically driven instantaneous spectral- 
model to the time-averaged spectrum may not be 'accept- 
able'QFor these reasons, judging from the best-fit residuals 
plot in Fig. 1 1 01 we believe that Model 4 provides a broadly 
acceptable fit to the 0.5-10 keV spectrum of the source. 

The new best-fit broad-band photon index is steeper 
than the one obtained from the best-fit Model 4 in the 3-10 
keV (Sect. [572. 3p . Interestingly, the new best-fit broad-band 
values for the remaining model parameters i.e. inclination 
angle, emissivity index, ionization parameter, and BH spin 
parameter, are entirely consistent (within the 90 per cent 
confidence limits) with those we have obtained from the 
best-fit Model 4 to the 3-10 keV spectrum. 

With respect to the BH spin-parameter, a, Fig. [11] 
shows the form of the A.\ 2 plot versus a, used in order to 
derive its uncertainties 8a, representing the 90 per cent con- 
fidence limit. Note here, a value of zero as well as the value 
of unity, are excluded at the 1 per cent significance level. 
We therefore conclude that relativistically blurred reflection 
from the innermost part of a slightly ionized disc, around 
a moderately spinning BH, may indeed be responsible for 
both the broad, asymmetric iron line and the 'soft-excess' 
emission that we detect from Fairall 9. 

After estimating the best-fit parameters of Model 
4 for the overall spectrum, we can compute the source 
X-ray flux and X-ray luminosity. In the 0.5-10 keV 
and 2-10 keV energy-ranges, the absorbed[unabsorbed] 
flux is (2.5[2.6]±0.2)xl0 -11 erg cm -2 s^ 1 and 
(1.5[1.5]±0.1) x 10 -11 erg cm~ 2 s _1 , respectively. The corre- 
sponding X-ray luminosities for the same energy-bands are 
(1.26[1.33]±0.03)xl0 44 ergs" 1 and (0.76[0.77]±0.05) x 10 44 
erg s _1 , respectively (the flux and luminosity errors 
correspond to one standard deviation in order to be 
easily comparable with previous observations). Finally, 



4 Obviously, it would be best to split the observation into smaller 
segments, e.g. three segments of 40 ks, and then fit the various 
models to the spectra of those segments. However, the source 
count-rate is not high enough in order to follow this approach 
and constrain the model parameters accurately. 
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Figure 10. Spectral fits in the 0.5-10 keV energy-range of Model 
4. Two narrow Fe lines — Fe Ka (fixed), Fe K/3 (fixed)-, a reflec- 
tion component (fixed), a redshifted power-law, and an ionized 
disc-reflection convolved with relativistic blurring model. In the 
ratio plot, the area between the horizontal lines, in the 0.5-2 keV 
energy-range, indicate the 8 per cent (dashed-lines) and 10 per 
cent (dottcd-lines) cross-calibration deviations between the pn 
and the MOS detectors. 



after applying the bolometric luminosity correction of 
iMarconi etalT (|2004l ) to the derived 2-10 keV X-ray 
luminosity, we get for Fairall 9 a bolometric luminosity of 
2.34xl0 45 erg s _1 . This value corresponds to an Eddington 
luminosity ratio of 0.07, comparable to the one derived by 
IWoo fe Urrvl (|2002T l (see Sect, [p but smaller by a factor of 
two than that derived bv lSchmoll etU] \200$ ). 



6 SUMMARY AND DISCUSSION 

We have analysed a long observation of Fairall 9, obtained 
with the XMM-Newton observatory during December 2009. 
EPIC light-curves have been produced and studied in several 
energy-bands, a thorough X-ray spectral analysis has been 
performed on both the RGS and the EPIC data-sets. 

The hardness-ratio analysis has revealed a 'softer- when- 
brighter' source-behaviour which is typical for radio-quiet 
type I Seyfert galaxies. The observed spectral variations do 
not have a large amplitude, In particular, the 5-10 keV light- 
curve (Fig. [2] lower panel) shows a small amplitude flux- 
increase of the order of 10 per cent and a small fractional 
variability amplitude, 2.2±2.2 per cent. This result implies 
that our best-fit model results for the 3-10 keV spectral re- 
gion should not be strongly affected by significant spectral 
variations which could make the study of the average spec- 
trum problematic. 

However, this is not the case for the lower energy-part of 
the spectrum between 0.5-2 keV (i.e. the 'soft-excess' spec- 
tral region) where the 0.5-1.5 keV light-curve (Fig. [2] top 
panel) exhibits an overall flux increase of 30 per cent and a 



10 D. Emmanoulopoulos et al. 



Table 1. Best-fit and fixed parameters for the X-ray spectral Models 2, 3 and 4 in the 3-10 keV energy-range 
and for Model 4 (broad-band) in the 0.5-10 keV energy-range. 



Number 


Model-parameters 


TV J ] „l o 

Model 2 


TV T„ J „1 Q 

Model 3 


^| ] „T ^ 

Model 4 


Model 4 (broad-band) 




Rcdshift, z 


0.047 


0.047 


0.047 


0.047 


1 


wabs: Column density, TVh (cm -2 ) 


3.20X10 2 " 


3.20X10 20 


3.20 xlO 20 


3.20X10 20 


2 


zGauss: Fc Ka line-energy (keV) 


6.42±0.01 


6.42 


6.42 


6.42 


3 


zGauss: Fc K/3 line-energy (keV) 


7.06 


7.06 


7.06 


7.06 


4 


zpowerlw: Photon index, T 


i oc+0.03 
± - oo -0.02 


1 97+°- 07 

-0.06 


1 84+ - 02 
1 -°^-o.oi 


2 01+ 01 


5a 


pexrav: Photon index, r pe x 


1.87±0.03 


1.87 


1.87 


1.87 


5b 


pexrav: Relative reflection, i? pex 


n O7+0.05 


0.87 


0.87 


0.87 


5c 


pexrav: Energy cut-off, E c (keV) 


200 


200 


200 


200 


5d 


pexrav: Elem. Z>2 abundance 


1 


1 


1 


1 


5e 


pexrav: Iron abundance 


0.8 


0.8 


0.8 


0.8 


5f 


pexrav: Inclination, <j> (degrees) 


40 


40 


40 


40 


6a 


laor: Line-energy (keV) 




6.75±0.18 






6b 


laor: Emissivity index, q 




o /ic+0-40 
2 - 46 -0.37 






6c 


laor: Inclination, 9 (degrees) 










6d 


laor: Inner radius (r g ) 




6 






6c 


laor: Outer radius (r g ) 




300 






7a 


kerrconv: Emissivity index, q 






2.071°;" 




7b 


kerrconv: Inclination, 6 (degrees) 






soil? 


64l^ 


7c 


kerrconv: Spin parameter, a 






43+0 07 


39+0-48 


7d 


kerrconv: Inner radius, (nsco) 






1 


1 


7e 

8a 


kerrconv: Outer radius, (risco) 
ref lionx: Iron abundance 






400 
0.8 


400 
0.8 


8b 


reflionx: Photon index, T 






Tied to param. 4 


Tied to param. 4 


8c 


reflionx: Ionization parameter, £ 
(erg cm s~ 1 ) 






IS.371^ 21 


15.82l^2g 


X 2 /d.o.f. 




600.29/538 


559.18/536 


578.88/535 


1052.63/749 



Values without uncertainties correspond to fixed fitting parameters. 



8 



6 



2 









99% 






1 95% 






I 90% 









0.2 0.4 0.6 0.8 1.0 



Spin parameter a 

Figure 11. A\ 2 plot versus the spin parameter a for the best- 
fit broad-band Model 4. The 90, 95 and 99 per cent confidence 
intervals for a are defined as the intersection points of the plot 
with the constant lines corresponding to Ax 2 of 2.706, 3.841 and 
6.635, respectively. 



fractional variability amplitude of 7.3±0.1 per cent. This in- 
dicates that any model fit to the overall broad-band average- 
spectrum will be affected by these variations. The 0.5-2 keV 
data have smaller uncertainties with respect to those in the 
2-10 keV, and thus the soft-band data drive the overall fit 
by having larger statistical weight. In fact, these differences 
in the amplitude variations between the soft and the hard 
X-ray bands can possibly account for the statistically low- 
quality fit of Model 4 to the overall spectrum of the source. 

The main results from our X-ray spectral analysis can 
be summarized as follows: 



(i) Our view to the central source in Fairall9 does not 
appear to be affected by significant absorption either from 
neutral or ionized material. 

(ii) Apart from a strong narrow iron emission-line and 
an absorption edge, we have also detected a broad emis- 
sion feature in the 6.1-7.7 keV energy-band, which can be 
well explained with relativistically blurred reflection from 
the innermost part of a moderately ionized disc, around a 
moderately spinning BH. 



Below, we discuss some implications of these two main 
results, within the context of current theoretical models of 
the X-ray emission in AGN. 
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6.1 On the absence of intrinsic absorption in 
Fairall 9 

Early observations by EXOSAT and GINGA showed 
that around 50 per cent of all type I Seyfert galaxies 
show strong evidence of warm-a bsorption (I Turner fc Pounds! 
1 19891 ; INandra fe Pounds! 1 1994 ). Moreover ASCA observa- 
tions sugges ted an even higher fraction of ar ound 50- 
70 per-cent (|Revnoldsl Il997l : iGeorge et al.lfl99sh . Finally, 
IBlustin et alJ l|2005l ). based on an XMM-Newton sample of 
23 type I Seyfert, found 17 objects having intrinsic ionized 
absorption, probably originating in outflows from the dusty 
torus. 

Nevertheless, there is a significant population of type 
I Seyfert galaxies that does not show any indication 
of strong X-ray war m-absorption features like: Aknl2 
JVaughan et all I2004T) Akn5 64 (iPapadakis et al.l I2007TI. 
Mrk 205 iReeves et all 120011). MA 335 ijGondoin et aD 
|2002|) . Mrk 4 7 8 dMarshall et all |2003| ) MCG-2-58-22 
(I Weaver et ail 1 19951 ), PKS 0558-504 JPapadakis et al 



20101) . and TON S 180 l| Turner et a~l l200ll: IVaughan et a 



2002) (but for this source see also IRozahska et alJ (2004) 



for very weak absorption features). 

Our main result, from the spectral analysis of the RGS 
data, is that there is no indication of significant X-ray warm- 
absorbing material along the line of sight to the central nu- 
cl eus of the Faira l l 9. Th is is in agreement with the results 
of lGondoin et all (|200ll ). Even the 99.73 upper confidence 
limit of N H , 3.72 xlO 20 cm" 2 (Fig. [5]), is smaller then the col- 
umn density of the warm-absorbing mat erial in most type I 
Seyfert galaxies (e.g. IBlustin et al 120051 ). This implies that 



our observations are not 'contaminated' by additional ma- 
terial between the central engine and us, making Fairall 9 
an ideal source for studying the physical properties of the 
immediate environment of the BH. 

The reason for the lack of warm-absorbing material in 
Fairall 9 (and the other sources listed above) remains still 
not fully understood. One possibility is that a significant 
column density of ionized gas exists but it is too highly ion- 
ized to show significant spectral features. We marginally de- 
tected weak emission-lines from ionized oxygen for which 
log 10 (17) = —1.1. This is a rather small value, and certainly 
not indicative of the presence of highly ionized material in 
the vicinity of Fairall 9's nuclear region. Another possibility 
is that the covering fraction of the absorbing material is less 
than unity. This could naturally arise in the case when the 
absorbing material is associated with a disc- wind which rises 
perpendicularly f rom the disc at a narrow range of radii (e.g. 
lElvis et aLlkoOOh . 

Both the significant outflow velocity of the emitting 
well its main source of ionization, as inferred by the 
OvuKa (f/r) ratio, 4.3I24 (i- e - photo-ionization) , suggest 
that the emitting gas in Fairall 9 could be probably part 
of a nuclear wind driven out by radiation pressure. This 
same evidence is f ound in numerous oth er Seyfert galax- 
ies, both absorbed jGuainazzi et al.ll2008h and unabsorbed 
(Bian chi et alJ 12009). It is commonly interpreted as a bi- 
conical outflow, launched from the innermost region of the 
AGN and spatially co-located with the optical ionization 
cones fre quently seen in Qui emission in ty pe II Seyfert 
galaxies jOgle et al.ll2000l ; lBianchi et al.ll2006l ). 



6.2 The 'bare' central nucleus of Fairall 9 

Our X-ray spectral results from the analysis of the pn and 
MOS data of Fairall 9, are in full agreement with the 'canon- 
ical' picture for the central region in Type I AGN, as it has 
been established the last 20 years. For example, the best- 
fit power-law continuum slope-estimate, coming from fitting 
Model 4 to the broad-band data in the 0.5-10 keV energy- 
range, is r = 2. 01^n m - This va l ue is s lightly softer than the 
one derived by IGondoin et all (|200ll ). 1.80t°; 2 3 usin S ear ~ 
lier XMM-Newton data, (but still in accorda nce within the 
errors ), and harder from the one reported bv lSchmoll et al.l 
(|2009l ). 2 .09=1=0.01 using Suzaku data. This is absolutely in 
accordance with the spectral variability behaviour of the 
source as a function of flux. During this XMM-Newton ob- 
servations the source flux in the 0.5-10 keV energy-range is 
aroun d 2.5xl0~ n erg cm -2 s" 1 (Sect. [573|) . IGondoin et all 
ll200ll) derived an e qual value for the flux, in contrast to 
ISchmoll et ail (|2009h who derived a higher flux of 4xl0~ n 
erg cm~ 2 s" 1 . Therefore, given the 'softer-when-brighter' 
source-behaviour, it is not surprising that the XMM-Newton 
spectral slope is slightly harder than that derived from the 
Suzaku data. 

Moreover, we have detected a narrow emission iron line 
at 6.42 keV (in the source's rest-frame). The best-fit value 
of the eqw for this line, is 103=1 =19 eV which i s som ehow 
lower from th e ones derived by IGondoin et alJ l|200ll ) and 
ISchmoll et aT] l|2009l ). 120 eV (no errors are quoted) and 
130±10 eV, respectively, but still within the quoted errors. 
Nowadays, narrow iron lines, at 6.4 keV (in the source's 
rest-frame), are detected almost universally in bright AGN 
X-ray spectra. The eqw of this line in our spectrum is very 
consistent to t hose o bserved in other sources. For example, 
iBianchi et al.l ([2009) , studying the spectra of 81 AGN, esti- 
mated an average eqw of the narrow Fe Ka lin e of 76±6 keV. 
Simila r res ults have also been reported by INandra et alJ 
l|2007h and lde La Calle Perez et~ai1 (|2010l ). 

This narrow iron line is usually attributed to reprocess- 
ing from cold Compton-thick matter which is most prob- 
ably associated with the putative, distant molecular-torus 
in Seyfert galaxies (e.g. iMatt et alJ Il99ll) . Based on our 
best-fit results of Model 3, we find a reflection fraction of 
R — 0.86 ± 0.17. From our estimate, the distant reflector 
should subtend an apparent solid angle of around 1.77T sr 
at the X-ray source. However, a more accurate statement 
about this and about the opening angle of the torus, can 
be made only by fitting to the data a physical model of the 
torus itself, rather than a slab (as we have done by using 
the XSPEC model pexrav) to account for the continuum 
reflection associated with the narrow ir on line. 

Following the best-fit linear mode l (IBianchi et al.H200 
of th e Iwasawa-Taniguchi (IT) effect (|lwasawa fe Taniguch 
1993), indicating the scenario of anticorrelation between the 
narrow Fe Ka emission-line and the 2-10 keV luminosity, we 
expect an eqw of 56±4 eV (for a source luminosity in 2-10 
keV energy-range of 0.77xl0 44 erg s _1 , Sect. I5,3"1l . In order 
to check the consistency between our direct estimates and 
the value predicted by the IT effect, we perform a z-test. 
The value of ^-statistic is -4. having a chance coinci- 



5 Note that the standard deviation of the estimated eqw of Fe 
Ka line is 10 eV 
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dence probability of only 1.17xl0 -5 . That means that with 
more than 99.99 per cent confidence we can reject the null 
hypothesis that the measured eqw of the Fe Ka line is com- 
parable to that predicted by the IT relation. 

The broad residuals of Model 2 around 6.1-7.7 keV 
(Fig. top panel), can be modeled very well by a relativisti- 
cally blurred disc-reflection model. A similar broad-emission 
spectral featu re has also been obs erved by Suzaku from the 
same source (|Schmoll et al.| [2009). Due to the lack of any 
significant warm-absorption in the nucleus of Fairall9, this 
picture should unveil directly the physical properties pre- 
vailing in the immediate environment of the BH. What is 
even more interesting is the fact that the same relativisti- 
cally blurred disc-reflection model can also account satis- 
factorily, under certain assumptions (see Sect. HO)) , for the 
'soft-excess' emission that we also detect from Fairall9. 

The emissivity index that we have derived from the 
best-fit Model 4 to the overa ll X-ray spectrum is q = 
1 - QS -o° l - ISchmoll et ail |2009D found an upper bound of 
5 but more string ent constraints have been derived by 
Patrick et all (|201lh . 2.7l°l, which agrees reasonably well 
with our findings. 

Our best-fit value for the BH spin parameter for 



Fairall9, a = 0.39Z 



is in accordance with the one de- 



rived by ISchmoll etal} l|2009I l. a = 0.60 ± 0.07, using the 
Suzaku data. The latter value comes from a broader-band 
spectral-fit in the 0.7-10 keV (X-ray Imaging Spectrometer) 
and 12-30 keV (Hard X-ray detector/PIN) whereas ours has 
been derived from a smaller spectral-range of 0.5-10 keV. 
Note that after ignoring t he spectral measurements below 2 
~" l|20Qg| ) got a value of a = Q.St^g. Also, 



keV. ISchmoll et al 



IPatrick et al.l ( 2011 ) , used the same Suzaku data of Fairall 9, 
together with the hard X-ray data from the BAT, on board 
the SWIFT satellite. By performing spectral-fits in the 0.5- 
100 keV energy-range (excluding 2.7-1.95 keV and 10-15 
keV) they obtained a value of a — 0.44^q'jj which is even 
more consistent with our estimate. 

Recent theories studying the effects of chaotic accre - 
tion episodes on the BH spin-parameter (|King et al.| [2008). 
predict that a should peak at around 0.1-0.2 for a BH 
mass of 10 8 Mq, depending on the vertical viscosity of 
the accretion disc. On the other hand, prolonged accretion 
should resu lt in much larger spin-p arameters, of the order 
of a > 0.9 l|Berti fe Volonteril 120081 ). Consequently, the BH 
spin-parameter of around 0.4, that we found for Fairall 9, is 
more consistent with models which favour a chaotic accre- 
tion scenario. 

Up to now, rather few studies have been conducted 
to specify the BH spin-parameter, a, for the type I 
Seyfert galaxies. A value of a = 0.9891^ ^ has been 
repo rted for MCG-6-30- 15 (|Brenneman fe Reynolds] 120061 .1 
and (|Miniutti et alj|2009h found for SWIFT J2 1 27.4+ 5654 a 
value of a = 0.6 ± 0.2. Moreover IZoghbi et all (|2010T ) found 
a high BH spin parameter value of a > 0.98 for 1H0707- 
495 and iGallo etai] d2010l) derived a value of 0.7±0.1 for 
Mrk 79. Finally. IPatrick et al.l (|201 lh suggested intermediate 
BH spin-parameters for Mrk 335 and NGC 7469 of 0.70t£;Ji 
and 0.69±0.09, respectively. 

This long-look XMM-Newton observation of Fairall 9 
has revealed that we are dealing with a 'clean' source since 
our X-ray spectrum is not affected by significant warm- 
absorption. Having this uncontaminated view of the central 



engine, a single relativistic disc-reflection model can explain 
the overall emission from 0.5 keV up to 10 keV. The exis- 
tence of a moderately spinning Kerr BH is strongly favoured 
from our data consolidating firmly that this source is in the 
intermediate spin-class of AGN. 
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